Constraints on the IR behavior of the gluon propagator in Yang-Mills theories 



A. Cucchieri 

Institute/ de Fisica de Sao Carlos, Universidade de Sao Paulo, 
Caixa Postal 369, 13560-970 Sao Carlos, SP, Brazil 



o 
o 

(N 

o 

Q 

o 

(N 



T. Mendes* 

DESY-Zeuthen, Platanenallee 6, 15738 Zeuthen, Germany 
(Dated: February 2, 2008) 

We present rigorous upper and lower bounds for the zero- momentum gluon propagator D(0) of 
Yang-Mills theories in terms of the average value of the gluon field. This allows us to perform a 
controlled extrapolation of lattice data to infinite volume, showing that the infrared limit of the 
Landau-gauge gluon propagator in SU(2) gauge theory is finite and nonzero in three and in four 
space-time dimensions. In the two-dimensional case we find D(0) = 0, in agreement with Ref. [lj]. 
We suggest an explanation for these results. We note that our discussion is general, although we 
only apply our analysis to pure gauge theory in Landau gauge. Simulations have been performed 
on the IBM supercomputer at the University of Sao Paulo. 
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I. INTRODUCTION 

Color confinement is a basic feature of hadron physics 
that still lacks a clear theoretical understanding. Among 
several explanations suggested in the literature (see [2j 
for a recent review) ,the so-called Landau-gauge Gribov- 
Zwanziger scenario [3, 4] relates gluon confinement to the 
infrared (IR) suppression of the gluon propagator D(p 2 ), 
whereas quark confinement is related to the IR enhance- 
ment of the ghost propagator G(p 2 ). This scenario is sup- 
ported by several studies using functional methods 0] . In 
particular, these studies @, 0, HI predict, for small mo- 
menta, a gluon propagator D(p 2 ) oc p 2 '""" 1 ' and a ghost 
propagator G{p 2 ) oc l/p 2 ( 1+aG \ The IR exponents od 
and ac should satisfy the relation ao — 2a>G + (4 — d)/2 , 
where d is the space-time dimension and should have 
a value in the interval [1/2, 1]. Clearly, if ao > 1 one has 
D(0) = 0, implying maximal violation of reflection 
itivity Q. In the four-dimensional case one finds 
ac ~ 0.59 and ao = 2ac- Similar power behaviors have 
also been obtained for the various vertex functions of 
SU(AQ Yang-Mills theories @, ©]. As a consequence, 
the running coupling constants from the ghost-gluon, 
three-gluon and four-gluon vertices are all finite at zero 
momentum, displaying a universal (qualitative) behavior 
Let us note that a key ingredient of these results is 
the non-renormalizability of the ghost-gluon vertex, i.e. 
Zi(p 2 ) = 1, which has been verified at the nonperturba- 
tive level [lj| using lattice Monte Carlo simulations. 

One should stress, however, that different IR behav- 
iors for the Landau gluon and ghost propagators have 
also been proposed in the literature. For example in 
Ref. [ll[ the authors find that -D(O) is finite and nonzero 
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and that ac ~ 0, with a gluon propagator character- 
ized by a dynamically generated mass. Similar results 
are obtained in [12]. On the other hand, in Refs. us- 
ing Ward-Slavnov- Taylor identities, the authors conclude 
that D(p 2 ) should be (probably very weakly) divergent at 
small momenta and that ac — 0. Recently, Chernodub 
and Zakharov [lj) obtained the relation 2ao + ao = 1 
for the 4d IR exponents of gluon and ghost propagators, 
by considering the contribution of these propagators to 
thermodynamic quantities of the system, such as pres- 
sure and energy density. This result, together with the 
previous relation between on and ac implies ao = 2/5 
and etc = 1/5, i.e. the ghost propagator blows up faster 
than p~ 2 at small momenta, while the gluon propaga- 
tor diverges as p~ 6 / 5 . Very recently, in Ref. |15|, it was 
shown that using the Gribov-Zwanziger approach one can 
also obtain a finite -D(O) gluon propagator and ac = 0. 
Finally, phenomenological tests [16| seem to favor a finite 
and nonzero D(0). 

Numerical studies using Monte Carlo simulations sug- 
gest that the gluon propagator is finite at zero momen- 
tumfl . IVn . [H, QjJ H0| and that the ghost propagator 
[E miHl] is enhanced when compared to the tree-level 
behavior p~ 2 . Moreover, in 2d and in 3d [J,[2(| the gluon 
propagator D{p 2 ) shows a maximum value for p of a few 
hundred MeV and decreases as p goes to 0. On the other 
hand, in 4d, even using lattices with a lattice side of 
about 10 fm, one does not see a gluon propagator de- 
creasing at small momenta [ijj]. It has been argued that 
an IR decreasing gluon propagator can be obtained nu- 
merically only when simulations are done on large enough 
lattice sizes |22j. However, from recent studies in 4d us- 
ing very large lattice volumes [H, [H, HH , one sees that 
D{p 2 ) either displays a plateau for momenta p < 100 
MeV or gets slightly suppressed at small momenta. Let 
us note that one of the main problems of the numerical 
studies of the gluon propagator is the lack of a simple 
way of extrapolating the data to infinite volume. 
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Here we discuss the behavior of the gluon propagator at 
zero momentum. We show that, instead of studying D(0) 
directly, it is more convenient to consider the quantity 



M(0) 



1 



dim 



(j . b 



(1) 



In a spin system this would be equivalent to studying the 
average absolute value of the components of the magne- 
tization instead of the susceptibility, which is of course 
a much noisier quantity. (Note that, by symmetry, the 
field components will average to zero if no absolute value 
is taken.) In order to relate M(0) to D(0) we derive 
in Section HT1 rigorous lower and upper bounds for -D(O), 
which are expressed in terms of M(0). Numerical data 
are obtained from extensive simulations in two, three and 
in four dimensions, for the pure SU(2) case, using very 
large lattices in the scaling region. We show in Section lllTl 
that using these bounds for D(0) and with present lattice 
sizes we have clear control over the extrapolation of the 
data to the infinite- volume limit. In the same section we 
suggest a possible explanation of the results obtained. Fi- 
nally, in Section HVl we present our conclusions. We note 
that our discussion concerning the bounds for the gluon 
propagator is general, although we only consider here the 
Landau-gauge propagator and pure SU(2) gauge theory. 
Note also that recent studies [2^, [26| have verified the 
analytic prediction that Landau-gauge gluon and ghost 
propagators in SU(2) and in SU(3) are rather similar. 
Thus, we expect that the analysis presented here should 
apply also to the SU(3) case. 



II. LOWER AND UPPER BOUNDS FOR D(0) 

As said in the Introduction, interesting lower and up- 
per bounds for the gluon propagator at zero momentum 
-D(O) can be obtained by considering the quantity M(0) 
defined in Eq. fT}, i.e. the average of the absolute value 
of the components of the gluon field at zero momentum. 
These components are given by 



(2) 



In Refs. "4] it was shown that in Landau and in Coulomb 
gauge the quantity M(0) should go to zero at least as 
fast as 1/N in the infinite- volume limit, where N is the 
number of lattice points per direction. This result is sim- 
ply a consequence of the positivity of the Faddeev-Popov 
matrix, i.e. it applies to gauge-fixed configurations that 
belong to the interior of the first Gribov region. 

In order to find the lower and upper bounds for D(0) 
lets us consider the inequality 
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where a; is a vector with m components Xi. This result 
simply says that the square of the average of an observ- 
able is smaller than or equal to the average of the square 
of this quantity and is equivalent to the inequality 
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Note that expression ([3]) becomes an equality when Xi = 
constant. We now apply ([3]) to the vector with m = 
d{N 2 - 1) components {\A b {Q)\). This yields 



(5) 



Then, we can apply the same inequality to the Monte 
Carlo estimate of the average value 



(6) 



where n is the number of configurations. In this case we 
obtain 



(I4(°)D 2 < (I4(°)i 2 

Thus, by recalling that 

V 



(7) 



(8) 



and using Eqs. ([5]) and (J7J) we find 

V(M{0)) 2 < D(0) . (9) 
At the same time we can write the inequality 



jU . b 

This implies 



;Ei4(°)i 2 ) ^ ({Ei4(°)i} 2 ) 



d(o) < vd{N 2 c - 1) (M(ay 



(10) 



(11) 



Thus, if M (0) goes to zero as V~ a we find that D(0) -> 0, 
< D(0) < +oo or D(0) — > +oo respectively if the 
exponent a is larger than, equal to or smaller than 1/2. 
Finally, let us note that the inequalities © and (fTTj) can 
be immediately extended to the case D(p 2 ) with p / 0. 



III. RESULTS 

We have considered several lattice volumes in 2d (at 
(3 = 10, up to a lattice volume V = 320 2 ) in 3d (at 
13 = 3, up to V = 320 3 ) and in 4d (at (3 = 2.2, up to 
V = 128 4 ). Details of the simulations will be presented 
elsewhere [13] ■ We set the lattice spacing a by considering 
the input value a 1 / 2 — 0.44 GeV, which is a typical value 
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for this quantity in the 4d SU(3) case. Note that the 
lattice volumes 320 2 at = 10, 320 3 at = 3 and 128 4 
at (3 = 2.2 correspond (respectively) to V « (70 fm) 2 , 
V « (85 fm) 3 and to V « (27 fm) 4 . Simulations in 2d 
have been done on a PC cluster at the IFSC-USP (with 
4 PIV 2.8GHz and 4 PIV 3.0GHz). Simulations in 3d and 
in 4d have been done in the 4.5 Tflops IBM supercom- 
puter at USP [28[. The total CPU-time was equivalent 
to about 5.7 days (in 3d) and 25.9 days (in 4d) on the 
whole machine. 

We start by considering the quantity (M(0)). We find 
(see Fig. □ and Table 0} that our data extrapolate very 
well to zero as 1/L 1 , with the values of I given in Table 
H Thus, in 3d and in 4d we have (M(0)) - l/V 1 / 2 , 
implying D(0) > 0. In particular, from our fits we obtain 
D(0) > (Bi/a 1 ) 2 . This gives D(0) > 0.4(1) (GeV~ 2 ) in 
3d and D(0) > 2.2(3) (GeV~ 2 ) in 4d, where we used 
a = 1.35687 GeV^ 1 in 3d and a = 1.066 GeV" 1 in 4d. 
As for the upper bound (|11|). using our fits (see again Fig. 
□ and Table© we have D(0) < d(N 2 - l)B u /a u , yielding 
D(0) < 4(1) (GeV" 2 ) in 3d and D(0) < 29(5) (GeV" 2 ) 
in 4d. On the other hand, in 2d both the lower and the 
upper bounds extrapolate to zero, implying D(0) = in 
agreement with Ref. Q. Let us note that our bounds 
in 3d and in 4d are in agreement with the data shown 
in Figs. 1 and 2 of Ref. [23]. [In the 3d case, c omp ared 
to the extrapolation reported in Fig. 1 of Ref. [23| , one 
should also include here a factor (3 — 3.0, i.e. 1.2(3) < 
D(0) < 12(3)]. Also note that in the three cases one 
finds B u « Bf and u 21. Indeed one can check that 
(M(0)) 2 < (M(0) 2 ), implying that the quantity M(0) is 
almost the same for all Monte Carlo configurations. More 
precisely, we verified for the three cases that ( M (0) 2 ) — 
(M(0)) 2 [i.e. the susceptibility of M(0)] goes to zero as 
~ l/V in the infinite- volume limit. 

In order to interpret these results, les us first note 
that, given a Gaussian random variable x with null mean 
value and standard deviation a, the random variable |x| 
has mean value (and standard deviation) proportional to 
a. In our case, this suggests that the average value of 
the gluon field at zero momentum A(0) [defined in Eq. 
@] should be zero with a standard deviation of the or- 
der of 1/L C , with c w I (see Table UJ. This is indeed 
the case in 2d, 3d and in 4d. [We find respectively 
c = 1.36(2), 1.47(3), 1.97(1) for the three cases.] In 3d 
and in 4d our results imply a oc 1/ \fV . This property is 
known as self-averaging [29| and is the behavior expected 
for extensive quantities in pure phases, away from phase 
boundaries. (In our case the magnetization is not exten- 
sive because we divide by the volume, but the result holds 
for the relative standard deviation.) More precisely, one 
talks of strong self-averaging when a oc 1/ L c and c = d/2, 
and of weak self-averaging when c < d/2. Thus, we find 
strong self-averaging for M(0) in 3d and in 4d and some 
kind of over self-averaging in 2d, with c > d/2. In sim- 
pler terms, the gluon propagator may be thought of as 
the susceptibility associated to the magnetization M(0) 
[or rather to the quantity defined by Eq. (TTJ) without the 
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FIG. 1: The square of the quantity o(M(0)) and the quantity 
a 2 d(N? - 1){M(0) 2 } (both in GeV -2 ) as a function of the 
inverse lattice side 1/L (GeV) for the 2d case (top), 3d case 
(center) and the 4d case (bottom). We also show the data 
for D(0)/V (also in GeV" 2 ) and the fit of the data using the 
parameters reported in Table [U 



absolute value, which has zero average]. In 3d and 4d 
the system has (finite) nonzero susceptibility, while for 
2d the susceptibility is zero. We do not have a simple 
explanation for this latter result. Here we can only argue 
that the 2d case is probably different since there are no 
propagating degrees of freedom. 

Note that our results in 3d and in 4d only imply that 
reflection positivity is not maximally violated. A clear 
violation of reflection positivity [27], [30( is still obtained 
in 2d, 3d and in 4d for the SU(2) and SU(3) cases. 
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TABLE I: Fits of a(M(0)>, a 2 (M(0) 2 } and D(0)/V respec- 
tively using the Ansatze Bi/L 1 , B u /L u and B/L k . Note that 
Bi, B u and B have mass dimensions respectively —l—l, —u—2 
and — k~ 2. Note also that in order to obtain Fig.[Tjone should 
multiply by d(N 2 — l) the data and the fit related to the fourth 
and fifth columns of the table. Most of the data used for the 
fits have a statistical error of the order of 2-3 %. For all fits 
we have x 2 /d.o.f. « 1. 

d Bi I B u u B k 

2 1.48(6) 1.367(8) 2.3(2) 2.72(1) 3.3(2) 2.73(1) 

3 1.0(1) 1.48(3) 1.0(3) 2.95(5) 1.5(3) 2.96(4) 

4 1.7(1) 1.99(2) 3.1(5) 3.99(4) 4.7(8) 3.99(4) 

IV. CONCLUSIONS 

We have shown that the Landau-gauge gluon propa- 
gator at zero momentum D(0) is finite and nonzero in 
3d and in 4d. At the same time, we find D(0) = in 
2d, in agreement with Ref. These results have been 
obtained by considering the inequalities in Eqs. ^ and 



(|llj) . i.e. by studying the "magnetization- like" quantity 
M(0) instead of the "susceptibility" D(0). This allows 
control of the extrapolation of the data to infinite vol- 
ume. Moreover, the quantity D(0)/V can be well fitted 
in this limit as a function of 1/L. Our results in 3d and 
in 4d can be explained as a manifestation of strong self- 
averaging. As mentioned above, a similar analysis may 
be applied to more general cases, and considering also 
nonzero momenta. 
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